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Abstract 

We propose a way of a dynamical solution of the strong CP in models with compact extra- 
dimensions. To this aim we consider a one dimensional toy model for QCD, which contains a 
vacuum angle and a strong CP like problem. We further consider a higher dimensional theory, 
which has a trivial vacuum structure and which reproduces the perturbative properties of the toy 
model in the low-energy limit. In the weak coupling regime, where our computations are valid, we 
show that the vacuum structure of the low-energy action is still trivial and the strong CP problem 
does not arise. Also, no axion-like particles are generated in this setup. 
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I. INTRODUCTION 

Quantum Chromodynamics has a complicated vacuum structure Due to the non- 
trivial group of mappings from the coordinate space to the gauge group tt^{SU{3)) = Z, 
there is an infinite number of perturbative vacua \n) which are not connected by smooth 
gauge transformations. Physical passage between these vacua is possible at an exponentially 
suppressed rate. One can further define physical vacua, that is to say, vacua between which 
transition is not possible. These are called ^-vacua and are expressed as discrete Fourier 
transform of the n-vacua, \6) = Xl^-oo '^^^ physics depends on which vacua we 

are in, so 6 becomes a new parameter of the theory. Equivalently this may be expressed 
by a new effective term in the Lagrangian of the form i6q{x), where q{x) is the topological 
charge density.^ The presence of this term implies that {q{x)) does not vanish and is an 
odd function of ^ ^]. The non- vanishing of {q{x)) leads in turn to the CP- violation and can 
be measured experimentally. In principle 6 can take any value in the interval [0,27r), but 
according to the measurement of the neutron dipole moment, it turns out to be extremely 
small 6 < 2 ■ 10~^°. The fact that this value is unnaturally small is known as the strong CP 
problem. 



Different models were proposed to explain the smallness of 6. Some scenarios pjj lead to 
a larger probability of living in a universe with small 6, but most of the models rely on the 
addition of some dynamical properties to 6. A possible way is to promote the parameter 6 
to a field, the axion, which will relax dynamically to zero [3| . This solution is however, in its 



minimal setup, ruled out by the non-detection of this axion particle [s]. More complicated 
models exist and can avoid these bounds, some of these also use extra-dimensions and allow 



the axion to propagate in the bulk j^j. 

Though both the CP problem and its axion solution have their origin in the symmetry 
properties of the system (existence of large gauge transformation for the first, and Peccei- 
Quinn symmetry for the second), one should be careful with the high energy extensions of the 
theory. In case where a high energy effect break some of the symmetries, the situation may 
change. This, for example, happens in the case of the global Peccei-Quinn symmetry. It is 
usually considered that it should be broken by higher order operators arising from quantum 
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gravity corrections [12|. The contributions from these operators scale hke a power of fa/Mpi, 
where fa is the axion couphng constant. Due to the astrophysical bounds from supernova 
coohng by hght particles, fa has to be large and the ratio fa/Mpi cannot be small enough, 
leading to corrections to the potential that spoils the solution of the strong CP problem. 
These problems are avoided in invisible axion models 1^ , where the Peccei-Quinn symmetry 
is a consequence of gauge invariance, and thus survives the quantum gravity effects. If the 
high energy theory does not violate the Peccei-Quinn symmetry, the solution to the CP 
problem is unaffected. On the other hand, if the high energy theory breaks the invariance 
against large gauge transformation, this may remove the CP problem by itself. This is the 
essence of the solution of the CP problem proposed in this article. 

Extra dimensions can provide a solution to the strong CP problem without light axion 
particles 0, Q, Q]. As we will see, the solution in this setup is also dynamical, but does 
not imply the existence of some light particle or any additional degree of freedom. The 
basic idea is that the group of mappings of a higher dimensional space to the gauge group 
{TT'^{SU{3)),d < 4) is trivial^ and 9 vacua cannot be defined. This type of solutions were 
already proposed some time ago with compact extra dimensions ^| or with infinite ones 
Q]. It was however pointed out in [8] that the solution might be more complicated than 
expected from the simple topological argument. The presence of a suitable extra dimension 
removes the degeneracy of the \n) vacua, thus preventing to define ^- vacua. This solution is 
only formal and the strong CP problem may reappear when we require that the observable 
effects of the extra dimensions are small. In this limit it can be expected that the usual 
QCD is recovered and the strong CP problem reappear in some way.^ To show that the 
strong CP problem actually disappear, we should not only get a trivial vacuum structure, 
but also show that the physical consequences of a non-vanishing 6 also do disappear, that 
is to say, the topological charge should vanish in average ((g) = 0). 

To our best knowledge, no complete setup where fields are explicitly localised on a brane 
and where the strong CP problem is solved have been constructed till now. Solving the strong 
CP problem in real QCD using extra dimensions is difficult. The full QCD is complicated and 
the localisation of non-Abelian theories on a brane is non-trivial [9|]. In particular the gauge 



^ Except Tr*{SU{3)) — Z2, which should not be an issue for our discussion. 
This is indeed the case in [sj when the extra-dimension is taken to be an orbifold. 
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theories in more than 4 dimensions are generically not renormahsable (while SU{N) gauge 
;heories might however be nonperturbatively renormahsable in 5 or even 6 dimensions, see 



23|). Thus, an eventual extra-dimensional extension of QCD might need an UV completion. 



that is to say will itself have to be thought as an effective theory comming as a low energy 
limit of some renormahsable theory (for instance string theory). However, certain important 
features the solution of the strong CP problem can be addressed in a simplified model. We 
will consider here the simplest possible model where a strong CP like problem occurs and can 
be solved by the addition of an extra dimension. Keeping this model as simple as possible 
will enable us to perform the relevant computations explicitly, retaining all the important 
topological properties. 

The simplest model which contains a QCD like 6 angle is the Abelian Higgs model in 1+1 
dimensions. This model has the same complicated vacuum structure as QCD (7r^(f/(l)) = Z) 
and its effective action also contains a ^-parameter. This model was already studied many 
times as a toy model for QCD. For instance the QCD U(l) problem was first solved there 

n 

[131 ]. This solution could then be mapped to QCD pj(. Note that the Abelian Higgs model 
was also successfully used as a toy model for the electroweak baryogenesis HQ. Our 
calculations could also be relevant to such issues. 

We then consider a higher dimensional theory, which reproduce the perturbative prop- 
erties of the Abelian Higgs model as a low-energy action. In this latter theory, the spa- 
tial dimensions have the topology of a sphere, so that the vacuum structure is trivial 
(7r^(t/(l)) = 0). At low energy, particles are locahsed on the equator (which will be the 
brane) and the extra dimension extends towards the poles. To simplify the computations we 
use in practice a geometrically simpler realization of this topology, and consider a pancake 
(the sphere is flattened to two disks). 

This article is structured as follows. We first review the relevant properties of the Abelian 
Higgs model in 1+1 dimensions in Section 2. In Section 3, we consider the Abelian Higgs 
model in 2+1 dimensions and localise the fields with the help of a warp factor. We show 
that the 1+1 dimensional theory is correctly recovered as low energy effective action and 
discuss the validity of the classical Kaluza-Klein decomposition in some detail. In Section 
4, we discuss the non-linear sector of the theory. The 2+1 dimensional equivalent of the 
1+1 dimensional vacuum structure is presented. We show that the degeneracy of the 1+1 
dimensional \n) vacua is lifted and this new structure is interpreted as a potential for the 
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topological charge density. Finally the evolution of the topological charge density is consid- 
ered. At least in the weak coupling regime where our calculations are valid, it relaxes very 
fast to zero, which solves the strong CP problem. We conclude in Section 6 and discuss the 
extension to more realistic models. 



II. TOY MODEL 



Toy models to address the strong CP problem in the framework of extra- dimensions 
were already considered. In 8|, 2-1-1 dimensional electrodynamics was localised on a one 
dimensional brane. When two spatial dimensions form a sphere, the topological charge 
density {q{x)) is shown to get the effective action of a harmonic oscillator with a small 
frequency. This kind of dynamics does not completely solve the 6 problem since {q{x)) is 
slowly oscillating in time. Furthermore, electrodynamics in 1-1-1 dimensions is not really 
adapted to our purposes. It a is trivial theory and contains a strong CP like problem only 
when charged particles are added. In [lo| the effective low dimensional action from [s^ was 
reconsidered with charged fermions and the strong CP problem was addressed but without 
reference to some extra-dimensional scenario. A simple model that contains 6'-vacua and a 
strong CP like problem is the Abelian Higgs model in 1-1-1 dimensions. This is the model 
that we will consider and before discussing its embedding into a higher dimensional space, 
we recall some of its properties [3]. 



A. Basic properties of the Abelian Higgs model in 1-|-1 dimensions 

The action for the Abelian Higgs model in 1-1-1 dimensions reads: 

S = I dxdt (^-\f^.F^'' - V{<P) + ^ \D^<P\'^ , (1) 

where = — ieA^j^, and = {H + {;)e"^ is a complex scalar field with the symmetry 
breaking potential 

n0) = ^(l0p-^^)^ (2) 

Note that we use tilde everywhere to denote 1-1-1 dimensional variables. We will now discuss 
the perturbative properties of this action. We want to identify the degrees of freedom, 
derive their spectrum and interactions. As we will latter compare the effective action of 
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the 2+1 dimensional Abelian Higgs model localised on the boundary of the disk to the 1+1 
dimensional Abelian Higgs model, we will directly consider the space dimension to be a 
circle of length 27rR. We always consider 1/i? to be much smaller that the particle masses. 
To extract the physical degrees of freedom in the simplest way, we use the unitary gauge 
(a = 0). We replace the coordinate x ^ R6 and decompose the fields in partial waves 



H 



\^ pine An 

* n 



(3) 



with (A")* = and (if")* = if"", so that the fields A^, H are real. Substituting the 
partial wave expansion ([3]) into the action ([1]) we get 



S 



n\2 



1 



2\v'^ + 



n 



n\2 



Xv 



SttR 



J2 H'^H'^H'^H- 



■n—m—k 



m.k 



GF- 



(4) 



We see that after symmetry breaking, the physical spectrum contains a scalar field (Higgs) 
with the mass m\j = 2Xv'^ and its dispersion relation is, as expected, 0^ = + 'n? /B^. 
The gauge field action reads 



GF 



1 ■ Pv'^ 

' \n\2 ^ " 



1 



n\2 



m 



2R 



/27tR 



^(A^A- - A-AT)h- 



—m—k 



m,k 



(5) 



The field Aq is not dynamical and can be integrated out using its equations of motion, 

.2 X _ 



n 



i?2 



R 



zV 



^27tR 



k m.k 

Inserting the latter equation in the gauge field action ([5]) and rescaling 



A'^ 



T^T^ A^ 

e V 



(7) 
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we are left with (up to the terms 0{{evR) ^)) 

n \_ ^ ' * m 

-— — - } A1A\K^h~'^~'^~^ + higher order interactions 

m,k 

This is the action of a scalar field with the mass = ev. Integrating out Aq leads to 
some complications, such as 0{{evR)~'^) corrections in the interactions terms and further 
higher order (in the fields) interactions. In the following we will mainly perform semi- 
classical calculations. To insure that they are valid, we will consider only small couplings. 
Consequently, corrections in ([8]) are small and not relevant to our discussion.'^ Note that the 
couplings in the Lagrangian ([T]) are dimensionfull. Weak coupling criterion means that 

A 

^ < 1 ^ > 1 (9) 

Of course the QCD coupling is large and this toy model is not realistic in this sense. However 
the strong CP problem does not rely on the fact that QCD is strongly coupled and we may 
very well use a simple toy model which enables simple perturbative calculations. 

As emphasised in the introduction, the Abelian Higgs model has almost identical non- 
perturbative properties as QCD. We list them in the following and, when it happens, point 
out where they differ from QCD. 



B. Vacuum structure 



In short, the ensemble of vacuum configurations are the ensemble of mappings from the 
space to the gauge group U{1). This ensemble can be divided into equivalence classes 
identifying vacua that can be related by a smooth gauge configuration. Mathematically this 
is represented by the homotopy group 7r^(f/(l)) = Z. This means that there is a discrete 
infinity of classes that cannot be related by a smooth gauge configuration. An element of 
the class n is for instance 

0("^=^exp(^), 4"^ = ^, ^ = 0. (10) 

^ They are not needed to compare with the effective action for locahsed 2+1 dimensional fields since we 
will make similar approximations there. 
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These equivalence classes can be distinguished by the Chern-Simons number: 

Ncs = ^j dxA,{x). (11) 

It takes the value n when applied on a vacuum state of the equivalence class n. The transition 
from one equivalence class to another can formally be achieved by a discontinuous gauge 
transformation. An example of the gauge transformation that changes the Chern-Simons 
number by n is 

^" = (^) • (12) 

C. Sphaleron 

Physically, the transition between two vacua needs to go through a set of non-vacuum 
configurations that form an energy barrier (see Fig. [1]). For instance, the set of static field 
configurations 

(j)'^'' = ve^ [cos(7rr) + i sin(7rr) tanh(m/^x sin(7rr))] , 

< - ^. (13) 

form a path that goes from vacuum n = at r = to vacuum n = 1 at r = 1, minimising 
the energy of the intermediate configurations. The configuration of maximal energy 

2 

Esph = ^mnv'^ (14) 

is reached at r = |, and is called the sphaleron. It is relevant for the high temperature 
behaviour of the theory. Thermal fiuctuations can reach the required energy Esph and the 
system may pass classically between vacua. Although all the previous and further discussed 
properties of the Abelian Higgs model exactly match QCD, we should point that, strictly 
speaking, the sphaleron does not exist in QCD. The corresponding configurations in QCD 
contain a free scale parameter and their potential energy can be made infinitely small in the 
limit of large configuration size. However, at large scale, kinetic energy along the path grows 
and the transition rate is still suppressed. This small difference will not be of importance 
for our discussion. 
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FIG. 1: Sketch of the vacuum energy as function of the winding number n 



D. Instanton 

At small or vanishing temperature, the system can also tunnel from one vacuum to 
another. In quantum field theory, tunnelling is represented by instantons, which are solutions 
of the classical equations of motion in Euclidean space-time. The set of configurations ([T3|) 
can serve as Ansatz to compute the Euclidean action, which gives the leading information 
to compute the tunnelling rate. 

The proper instanton in this model, i.e. the field configuration that minimises the Eu- 
chdean action and describes the tunnelling between the states |0) and \n) is the Nielsen- 



Olesen vortex 



2l| with winding number 



e 



Q = ANcs = ^J e.^F^^'^d'x = n. (15) 

Parametrising the Euclidean space-time in polar coordinates (r, 6) the field configuration 
reads: 

(f)(r,e) = e^"7(r), (16) 

A\r,e) = e'^PA{r), (17) 

where r is the unit vector r = (cos^, sin 6') and e'^^ the completely antisymmetric tensor 
with e^^ = 1. The functions A and / have to satisfy the following limits: 

'■^0 \n\ 



f{r) — > cr 
fir) ™ V 



A{r) ^ 0, (18) 

A[r) --. 

er 



Introducing dimensionless variables 



the classical action of the instanton is 



' A^^A, f^^f, (19) 



' + fir)' + f{rf Uir) - ' 



z / 1 \ 2' 



+ ^ rW-^ =^^V5(/^)' (20) 



2 V' 

where /i^ = \/(? and the function -B(/i) is of order one, and depends weakly on its argument. 
E. 6'-vacua and topological charge 

We shall consider the non-local gauge transformation that changes the Chern-Simons 
number by one. It commutes with the Hamiltonian and is unitary. The operator f/^ can 
therefore be diagonalised simultaneously with the Hamiltonian and must have eigenvalues 
of modulus 1. The 6'- vacua are defined as a superpositions of the \n) states which are at the 
same time an eigenvector of with eigenvalue e*^: 

|^) = ^e'"V)- (21) 

n 

A physical transition between the ^-vacua is not possible. They form different sectors with 
different physical properties and ^ is a parameter of the theory. The topological charge 
density q = e^i.F'^'^ in a 6'- vacuum does not vanish in average 

(g) = SnKe-^-' sin 6, (22) 

with Sci the instanton action and the factor K ~ rriH takes quantum corrections into account. 



III. 2+1 DIMENSIONAL ABELIAN HIGGS MODEL ON A DISK 

We will now recreate the 1+1 dimensional Abelian Higgs model as the low energy effective 
theory for a 2+1 dimensional model. The perturbative properties of our low energy theory 
has to resemble closely the 1+1 dimensional Abelian Higgs model, but hopefully have a 
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different vacuum structure. We expect that this is possible since, as 7r^(f/(l)) = 0, no ^ 
vacua exists in the original 2+1 dimensional theory. 

As discussed in the introduction, we will suppose that two spatial dimensions are the 
surface of a pancake, where the low energy fields live on its boundary and the extra-dimension 
extends from the boundary towards the centre. This pancake can be adroitly sliced into two 
disks. By symmetry we can expand the fields into functions that are odd or even with 
respect to pancake flipping. In the following we will only consider one disk where the fields 
can have either Neumann or Dirichlet boundary conditions on its boundary. 

In 2+1 dimensions, we consider a complex Higgs field and an Abelian gauge field. The 
space time is a disk, which we parametrise with polar coordinates (t,6,r). We will also 

n n 

introduce a warp factor A(r) to localise the fields [9|, [191]. Note that we do not want to 
consider gravity here and the warp factor will be thought as coming from the coupling to 
some external classical field, which multiplies the whole action. 



S= drde dt^A{r) 



(23) 



with the Higgs potential V^(0) 
disk is 



and Dm = — ieAu- The metric of the 



M- 



V 



-V 



and the volume element ^/g = r. The action can be rewritten more explicitly as 



S = J drde dtA |i 



+ 



In the following, we will use 



Or 



er 



r 



(24) 



ID, 



A 



— r- 



V ) 



A(r) 



(25) 



The exact form of A is merely chosen to simplify the algebra. It is only needed that 
J drA{r) < oo and that A(r) decreases from the brane (r = R) towards the extra dimension. 
In the following we will consider the following relation between the parameters: M ^ rriH ~ 
mw 3> 1/-R. 
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A. Linearised theory and spectrum 



We decompose the Higgs field in polar coordinates as = (f + H)e^'^, and retain the 
quadratic part of the total (gauge field + Higgs) Lagrangian ; 



^0 



'-A(r) 



-F^s + rFl. - -Ff 



Or ^ er 



ridoa — eAof' — ridrcr — eArY (decr — eAg 

r 



-rA—H^ H 



r 



(26) 



with m'jj = 2Af ^. The linearised equations of motion read (with the rescaling a — ea) 

-- 0, (27) 



doidocr - Ao) + \drrA{dra - A^) + ^de^dea - Aq 
r A r 



— 9,rA(9oA. - drAo) + ^de{d^Ae - dgAo) + ml,{A^ - d^a) = 0, 
rA 

r A 

-9o(9oA0 - deAo) + —dr — {drAe - deAr) - m\^{Ae - decr) = 0, 
-do{doAr - drAo) - ^deidrAg - dgAr) - m1y{Ar - drcr) = 0, 



rA 



drArdr + ^dl -dl-m]j]H = 



(28) 
(29) 
(30) 
(31) 



Note that the phase a and the gauge fields Am decouple from the physical Higgs H. These 
two sectors are analysed separately. 

a. Higgs sector: After Fourier transform 



the Higgs field /i" satisfies the following equation of motion: 



1 



n 



—drArdr - — -mir + uj^] K'ir) = 0. 



(32) 



(33) 



The explicit solution of the equation fl33l) which satisfy the boundary conditions /i"'(r) oc r^' 
at r — > 0, is given in terms of the Laguerre function: 



/i"(r) oc e-"(*^+^«)(^)"^'%^_ _i(2rl]H), (34) 



with VIh = a/ + ni\j — uS^. We should further impose that either h!{r) (symmetric mode) 
or hir) (antisymmetric mode) vanishes at r = i?. For n = 0, the lowest mode is symmetric 



/ig = const and has the energy cuq = rxiH- 
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For n 7^ 0, there is a low energy symmetric mode that satisfy, in the hmit M S> mn 



the usual dispersion relation uJq = y mj^ + ^ + (9(e~^^^^) (see Appendix so that we 
have a Higgs particle on the brane, with mass mn and, comparing to (jl]), correct dispersion 
relation. 

The next mode is antisymmetric and has uji > M. Further modes have higher energy 
and build a discrete spectrum (labeled with k) which satisfy ujk+i — t^/t = ^ at large k (see 
Appendix [XI). 

b. Gauge sector: After Fourier transform 

A^it,e,r)= /■^e-^-*^e^"%(r), (35) 

n 

the equations of motion for a, Oq , , a" are 



1 tTl 

iuji-i0a - tto) + —dr^ridrcr - <) + -^{ina - ag) = 0, (36) 



-^dr (Ar(z0< + dra^)) + ^(Ti0ae + ^'«o) + (^o + ^0(^Ww = 0, (37) 

r A 

0^a^ + 71000 + '^'^r— {drag - ina"^) - (a^ - ina)mly = 0, (38) 

2 

XTt Til 

0^a^ - i0(9ra^ - ^SrOg - -^a" - - 9,.cr) = 0, (39) 

with = e^f^ the W mass. For the continuity of the fields we require that Og , a" vanish 
at the origin for all n, except for which might be constant. Note that these four equations 
are not independent in the general case. A possible resolution is to work in the unitary 
gauge (cr = 0). We extract Og from the equation (136|) and replace it in the equations (l38l 
[39]): 

-^dr-dr -u' + % + ml\a-g- ^< = 0, (40) 

-a,— A - ^ — + j < + — = 0. (41) 

The equation (1371) for Oq is not independent and can be dropped. The analysis is shghtly 
subtle and we have to treat all possible cases separately. 

Consider first = 0, the equations fj^ W\\\ decouple and read (for r < R) 

/ 9 Mr — 1 \ 

[w' - m^) a^(r) + J (a^)'(r) + (a^)"(r) = 0, (42) 

-ml. + <ir) + (^^^) «)'(r) (43) 

+ «)"(r) - 2M5{r - R)a^^ = 0. (44) 
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The presence of the (5-function in the equation (jUj) changes the boundary condition on the 
brane for the symmetric mode to (a")'(i? — e) = —Ma^{R) = — (a")'(-R + e), which makes 
the lowest energy mode for a" heavy {u > M). For a^, we have the general solution 

^ g-KAf+n^)^2^2 (2rnH.) , (45) 

with Vtw = a/ + rri^ — uS^. The lowest energy mode satisfy (ag)'(i?) = and is light, 
cUq = miy + 0(6^^*"^^) and the next mode has energy of the order M. 

We shall now consider the case n 7^ 0. Equations (HOl HTl) are coupled and should in 
general be solved numerically. However if we notice that the mode we are looking for lives 
on the brane and therefore should have 

F,e = -m<(r) + 9,a^(r) = 0, (46) 

the equations (HOl HTj) reduce to (for r < R) 

[-mi^ + - ^) «^(^) + (2^^ + ^) («^)'(0 + {<)"{r) = 0. (47) 

This is similar to the equation for the Higgs mode, with mn replaced by mw- The modes 
are therefore similar to the Higgs ones, 

oc e-''(^+^'^V"L2%^_^_i {2rnw) (48) 

and a" = —jidrag{r) have to be continuous at the poles and one the brane. The lowest 
energy mode satisfy (ag)'(-R) = a'^{R) = and has energy Uq = ■\J + 75- + C(e~^*^-^). 
As for the Higgs, the next mode has energy ui > M. Like for n = 0, no other low energy 
mode is found in the equations (HOl l4T]) . 

To summarise, up to corrections (9(e~^^^^) <^ 1 the low energy spectrum and the disper- 
sion relations exactly match the usual 1+1 dimensional Abelian Higgs model. 

B. Effective action 

We can now built the effective 1+1 dimensional action for the Kaluza-Klein (KK) modes. 
To simplify the algebra, we will remove all the (9(e~^^^^) suppressed terms. We work in the 
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unitary gauge and use the following decomposition for the fields: 



-.0000 



V + 



e^-'hUrmit), 



A, 



An 



k=0 n=— 00 
00 00 

in9 „n 



fc=0 n=— c 
^ 00 00 

k=0 n=—c 
^ 00 00 

^EE 



(49) 



And n 



<Ar)AUt), 



And n 



aUr)Al{t), 



* k=0 n=— 00 

with k labelling different KK modes, and n labelling the angular momentum. Note that 
we have to impose that the fields are real, that is to say (A")* = y4~", (if")* = if ^" and 
{al,y = alP, m* = h-\ 

To get the 1+1 dimensional action for the whole tower of the KK modes we substitute 
the KK expansion (H9| [50]) in the Lagrangian (l23!l . and integrate over the extra-dimension 
coordinate r. In the following, we will consider only the = mode in the KK expansion fH9|) 
to get the low energy effective theory. We will also neglect the back-reaction of heavy Kaluza- 
Klein modes on the low energy action, supposing that these corrections are suppressed. This 
assumption will be checked in the next section. 

c. Gauge kinetic term Combining the relations (!36l [50]) we get 



W n 



The expansion (H9|) is inserted in the gauge kinetic term 

1 



GF 



dr dO dtA- 



lp2 , p2 _ lp2 



Or 



dr 



(50) 



(51) 



We use the formulas (l50l |46|1 , separate the t and r dependencies and get after some straight- 
forward algebra (keeping only the = mode) 



Sgf = I j dtY^]^{\A-\ 



m 



2 \2 

w) 



\A 



n\2 



with 



dr — 
r 



<\--,\W{r)[ 



(52) 



(53) 



the dot means derivative with respect to time and the prime derivative with respect to r. We 
fix the normalisation of such that i = 1 which, with the dispersion relation 0^ = + 
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leads to 



Sgf= I dtJ2U\A-\ 



(nVi?2)2 



n|2 



d. Higgs kinetic term We shall now consider the Higgs kinetic part of the action, 

r\Do(j)\^-r\Dr(j)\^--\Dg<f)\' 



Sd = J dr dO dtA- 
Again, we insert the KK expansion (149|) . use the relations (150| H6ll and get 



S 



D 



n 



m'w\Ayia + -\Hy / dr|/i'^(r)|' - —-\Hyih + 5 



lint 



with 



A 

r 



2 •\2 

Vy/ l„n|2 



dr/\r[W)\r)\' + ^^\h-{r)\' 



and we normalise /i" such that 

j rfr|/i"(r)|2 = 1. 
The interactions between gauge field and Higgs read 



\ ^ n,k,m V 2vr J 



with 



2 



, rA 



nm 



R^0n R^0m " " nm 
e. Higgs potential The remaining Higgs potential is 



which is expanded to 



V 



dt 



dr dO dt 



^ 



-rA 



Xv 



J2 H''H"'H-''~"'Id 



n,m 



n,m,k J 



with 



ArK'K^h- 



' I 

1 



n urn uk u —n—m—k 



ArK^K^h^h 
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/. Low energy theory In Appendix [B] we compute that 

,2 



L = -1 



and to 0{{MR)-\ {mHR)~^), we have 



n 



4 = 1 



M J _ M J _ M r_M 



which gives us 



5*, 



and up to 0((Mi?)-\ (m^i?)'^), 



2 



1^0 



I Qint 



nint 



V2^ 



j^n j^m j^k j_j—n—m—k j^m j_j^k m—n—k 



(66) 



(67) 

(68) 
(69) 



(70) 



Comparing to (jH [8]), we see that up to at least 0{{MR) ^, {muR) ^), we get the same 
interactions as in the 1+1 dimensional Abelian Higgs with the identifications 



= e^M, = vyM, X = AM. 



(71) 



Note that these identifications are just what we would expect from dimensional analysis. 
We conclude that the theory (123!) successfully reproduces the 1+1 dimensional Higgs model 
on the boundary of the Disk as its low energy limit. 



C. Validity of the Kaluza-Klein expansion 

Our computations are purely classical, and up to now we haven't taken into account the 
corrections to the effective action coming from interactions with heavy modes. This goes 
beyond the scope of this article and we will only describe here a region of the coupling 
parameter space where our semi-classical computations can be trusted. 

The first problem that occurs in most warped scenarios, is that the couplings become 
large far from the brane. For instance, this can be seen if we rescale the gauge field like 
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cAm — > Am- The gauge kinetic term becomes ^FmnF^^ and we see that the effective 
charge is As A(r) becomes very small far from the brane, the effective charge becomes 
very large. This leads to large quantum corrections and large effective couplings between 
modes which probe the region with small A. This problem first concerns very energetic 
modes that live in the bulk and can probe the region with small A but having strongly 
coupled heavy modes may also lead to large corrections to the low energy interactions. 

g. Behaviour of the theory at very high energies We will compute here the effective 
interactions between heavy modes and show under which conditions interactions between 
light and heavy modes do not affect the low energy theory. We only show explicit calculations 
for the Higgs modes. The gauge fields case is similar. 

If we compute the effective theory for all the KK modes we get the terms in the action 
of the form (EH ESD 

yn,m,p,— n—m—p 
fci,fc2,fc3,fc4 Trm ttP TT-n-m-p 

SnR -"fci-"fc2-"A:3-"fc4 
The quartic coupling between KK modes is given by 

KmZT = ^l^^{tl K (0 j dr. (72) 

For KK modes with very high energy 0^ ^ M, this integral can be approximated as (see 
Appendix O) 

Tni,n2,n3,n4 ^ 2M_R V ' (7'\\ 

It can be shown numerically that this formula fits well for all ki > 1. We shall now study 
what are the consequences on these large couplings. 

h. Corrections to the nonrenormalisahle couplings The simplest correction to the ef- 
fective action is the appearance of a nonrenormalisahle Agiff coupling, which comes at tree 
level from the following diagram. 




The contribution from heavy modes will be of the order 



^e-E^MM^- (74) 
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For heavy KK modes, / > 1, we can use equations fl73l) and the relation t^fc ~ f|, which 
leads to 

The sum converges fast but the factor e^^^ is very large. If we want our classical calculations 
to be valid, we have to impose a strong restriction on the 2+1 dimensional couplings A,e^, 
so that Xq and other corrections are small, that is 

A, < Me-2^«. (76) 

This assumption will be considered in the following; it enables to perform semi-classical 
computations. Although the assumption fl76l) seems very restrictive, there is no obvious 
reason that the further results are only valid under this constraint. We suppose that they 
are also valid under more general assumptions, but we will not prove it here. The extension 



to larger couplings requires lattice computations [22 1. 



IV. NONLINEAR SOLUTIONS 

After discussing the perturbative spectrum of the theory we study solutions of the nonlin- 
ear equations of motion. In the usual 2+1 dimensional Abelian Higgs model — not localised 
on a brane — neither sphaleron nor instanton exists. There only exists solitons — Nielsen- 
Olesen vortices in two spatial dimensions. The vortices possess the property that, if we 
consider a closed curve C surrounding the vortex centre, the following relation is satisfied 

|-/A-dl = n, (77) 

where n is an integer called topological charge of the vortex. 

nn 

Processes discussed in [8|, |20[ can happen: starting from the three dimensional vacuum, 
a vortex-antivortex pair is created near the south pole. The vortex can move across the 
equator up to the north pole. The bosonic configuration of the final state looks on the 
equator (brane) like the n vacua. Indeed using the relation flTTl) . the Chern-Simons number 
is 

Ncs = 7:~ dxAi = n. 



27r „ 

From the point of view of someone living on the brane, a topological transition (indistin- 
guishable from an instanton transition) occurred and the vacua changed (from |0) to |1)). 
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The difference with the pure two dimensional theory is that the end point is not a vacuum 
state. It contains two vortices, which represent a large energy cost for the system. 

We shall now see what happens when the warp factor is introduced and the fields are 
localised. The main effect of the warp factor is to damp the action of field configurations 
localised far from the brane. We therefore expect that, at the classical level, the vortex is 
light close to the pole and heavier when crossing the brane. 

The 1+1 dimensional \n) vacua with n are not vacua any more in the higher di- 
mensional theory. They become resonances with the energy of two vortices with topological 
charge n and some suppressed decay rate. 



A. Vortex at the centre of the Disk 



With the warp factor, the vortex is repulsed from the brane and a stable solution of the 
equations of motion can only be found at the centre of the disk. In this section we study 
the energy of the vortex at the centre of the disk and find its dependency on the parameters 
of the model and on its topological charge n. The vortex has the same properties as the 
1+1 dimensional instanton that is to say (f){r = 0) = 0, rotational symmetry and winding 
number Q = n. It is also parametrised by equations (fT6l[T71) but the functions /(r), A{r) are 
deformed by the presence of the warp factor. The energy of the A{r), f{r) field configuration 
is 



E 



2 2 



dr 



A'^ + A 



rf' + ^{n-Ar + r^{f~f^-r 
r 2 



(78) 



with 



^ and their equation of motion read 



- ^drA— - fin - A) 
A r ^ ' 

2, .2/ 



--drArf + h{n - AY + rV(r - /i"')/ 



0, 
0. 



(79) 



In the following, we will solve numerically the equations (!79|) with boundary conditions (fT8|) 
as in ITj. There are tree dimensionless parameters, one for the gauge-Higgs fields /i, one 
for the warp factor M = and the size of the Disk R = \/\vR. The dependencies of the 
energy with respect of M, R and the topological charge n will be studied. For the sake of 
simplicity we will set = 1. 
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The dependency on the topological charge is a parabola E (x for fixed M, R. An 
example of this is shown in Fig. [2l We checked that for M = 0..10, R = 1..6, n = 0..10, the 
discrepancy between the fit oc and the points is much less that 0.1% 

The dependence on the radius size Ris E oc exp{—2MR), with some deviation from this 
behaviour at small values of R, In Fig. [3l we plot the value of 

_ E 

~ 7rt;2/i%'exp(-2M/?)' ^ 

We will now extract the M dependency. We use of the previous dependencies and plot 
1] for M = 0..10. In Fig. HI r/ is fitted with values of n = 0..9, and the large R limit is used. 
For large M the dependency is again a parabola E oc M^, and deviation from this occurs at 
small value of M. 

To sumarise, within the relevant region (M ^ 1 ^ 1/^) we have the following vortex 
energy 

E ~ 47ry (81) 
Coming back to dimensionfull parameters, this takes the form 

E ^ 47r^n^e-2^^^. (82) 
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FIG. 3: Parameter r] as function of the disk radius R, Other parameters are n = 1, M = 4 




FIG. 4: Parameter r] as function of the brane scale M, fitted with 4.013M^. Discrepancies from 
points to the fits only occur for very small M 

B. Vortex sitting on the brane: Sphaleron 

The other case of interest is when a vortex of topological charge one is located on the 
brane, that is to say (f){R, 6') = at some 6 = 6; this is the 2+1 dimensional analog of the 1+1 
dimensional sphaleron. The presence of the brane breaks the rotational symmetry of the 
vortex. The exact solution is found by integrating numerically the 2 dimensional equations 



19.6 
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FIG. 5: Energy of the different topological solutions as function of the brane mass-scale (for 
= !)• The energy is scaled so that the 1+1 dimensional sphaleron has energy 1. The energy 
of the vortex sitting on the brane meets the 1+1 dimensional sphaleron energy for large M. The 
accuracy of the results is estimated to be better than 0.1% 



arising from the variation of the energy functionaL Details of the numerical procedure are 
given in Appendix [Dl 

The results of the numerical simulation are shown in Fig. [51 We only consider the case 
muR ^ 1. As long as M > 1, the 2+1 dimensional sphaleron energy exactly matches the 
1 + 1 dimensional one. Note that if the KK mass is larger than the sphaleron negative mode 
oj- ~ 0.8m//, the properties of the sphaleron are expected to be part of the low energy 
theory and match the 1+1 dimensional result (JT4l) . 



C. Movement of the global mode in the potential - solution to the strong CP 
problem 

The topological charge can be seen as a field [2I, in the 1+1 dimensional Abelian Higgs 
model it has a periodic potential (see Fig. [1]). The addition of a compact extra-dimension 
breaks this periodicity (see Fig. [6l). 

We shall compare the different scales in our model. As usual we consider the inequality 
M ^ uiH ~ Tnw ^ 1/R together with the 1+1 dimensional weak coupling requirement 
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FIG. 6: Sketch of the minimal energy for a configuration with winding number n. The value in 
local minima correspond to the energy of two vortices at the centre of the disk and the maxima 
contains the contribution of the sphaleron energy. 



V ^ 1 and the 2 + 1 one (1761) . The sphaleron mass 

^ ^ ^2 » 1 (83) 

is much larger that the particle mass in the weak coupling regime and much larger than the 
vortex mass,^ 

'^^p'^ ^ ^ 1. (84) 

^vortex M^^C 

At first sight, the energy of the vortex E^{n), and therefore the energy En = 2E^{n) of the 
state \n) seems very small. However if we consider that the couplings (electric charge here) 
have to be tiny to be in the regime where our classical approximations are valid, we have^ 

!!!=Sfi = _J!£L^ » » 1. (85) 

The latter relation shows that vortex pairs can decay to particles and the system relaxes 
to the n = vacuum. In this case the strong CP problem is solved. Indeed the 1+1 



^ In fact, the full set of requirements for the parameters is the following <C e~2Mfl ^ ^^h. ^ \ The 
first inequality is the weak coupling of the 2+1 dimensional theory, the second is ((M]) and the third says 
that the KK modes are heavier, then the localised ones. 

^ In fact, even stronger inequality is true, mvortex ^ M. 



24 



dimensional instanton transitions which make the average of the topological density ((g) 7^ 0) 
non-vanishing are not anymore instanton transitions and can occur only if a large energy is 
available in the system. At low energy, for example for the measurement of the neutron dipole 
moment, a topological transition cannot occur, the topological density vanishes ((g) = 
since already (|g|) = 0) therefore 9 is effectively zero. 

In the strong coupling regime, if the mass of the vortex happens to be too small the states 
1 77.) will be almost degenerate and we expect that the strong CP problem reappear at some 
point. Indeed when the \n) are almost degenerate, the real vacuum will be a supperposition 
of several \n). Several different supperposition of these \n) states with different phases and 
slightly higher energy will exist and we expect to have some sort of discrete strong CP 
problem there. The precise analysis of these states is left for future work. 

V. CONCLUSION 

Extra- dimensions might give an appealing solution to the strong CP problem. The so- 
lution relies on a global field dynamic and does not imply the existence of a new particle. 
Unlike for the axion solution, it cannot be ruled out by astrophysical observations. How- 
ever, the processes with change of the topological number (and, at the same time, fermion 
number) proceed differently in the proposed setup. 

In this paper the 2+1 dimensional toy model was studied, which at low energy reduces to 
the 1+1 Abehan Higgs model, often used to model some QCD properties. If the inequalities 
'f^sph ^ 'fnvortex ^ mj^, m^y are satisfied, the 1+1 dimensional \n) vacua become, for n 7^ 0, 
resonances in the 2+1 dimensional theory. These resonances have a high energy and decay 
reasonably fast to particles. The system relaxes to the |?t, = 0) vacuum, which is the true 
vacuum of the 2+1 dimensional theory. The vacuum structure is trivial and the strong CP 
problem is solved in this case. We showed that these inequalities ( l83l 18^ 1851) are realised 
in the weak coupling regime fl76l) . We suppose that these relations will remain true in the 
strong coupling regime. This seems reasonable for the relations (l83l [8^ while for (!85l) . it 
seems to fail, since we used the weak coupling assumption. However in the strong coupling 
case, the quantum corrections are large and, unlike the classical action, are not suppressed 
by the warp factor at the vortex centre. Indeed quantum corrections are expected to be of 
order h independently of the normalisation of the action. It would be interesting to study the 
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strong coupling regime with lattice simulations to check whether this large extra-dimension 
solution to the strong CP problem also works at large couphngs. 

We choose in this paper a very simple geometry, which allows for simple analysis but 
which has the annoying property that the size of the universe is just 27r larger than the 
extra-dimension. This restricts the semi-classical calculations to very weak couplings. A way 
to avoid strong coupling problems would be to consider a more complicated geometry, which 
allows for unrelated universe and extra-dimension size. We could then have a very large or 
infinite space and a small enough extra-dimension. The set-up will be more complicated 
and may not enable a simple analysis of the nonlinear solutions. However as we have shown 
here, it seems that the localisation of the fields on a brane can be done independently of the 
bulk geometry. It is indeed an interesting result on its own that an Abelian Higgs can be 
localised on the boundary of a disk. 

Finally, we may speculate how the obtained results transfer to the real world QCD. We 
expect that the degeneracy of the n-vacua will be lifted by the extra-dimension and the 
new \n) states (or resonances) will have an energy monotonically increasing with n. The 
tunnelling between the n and n — 1 will still be suppressed by the instanton (or, strictly 
speaking, bounce) action. At very high energy, these bounce transitions are possible and 
will look like normal instanton transitions. If we can arrange the geometry so that the 
states (or resonances) |n 0) have sufficiently large energy, bounces will not contribute in 
a measurement of the neutron dipole moment. The topological charge density will be zero 
and there will not be any observable effect of a ^-angle, so that the strong CP problem will 
be solved. 
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APPENDIX A: SPECTRUM OF KALUZA-KLEIN MODES 



i. Higgs dispersion relation We derive the dispersion relation for the hght mode of the 
Higgs. To this aim, we have to impose the boundary condition on the brane, h'{R) = 0. 
Suppose that the low energy mode satisfy u ~ rrin <^ M. We can then use the asymptotic 
expansion for large argument of Laguerre function : 



L^iz) ™ sin(z/7r)r(A + z/ + l)z-^-^-^e' + ' ^ ' + 0{z-')j . (Al) 

In this limit, equation (IMI) becomes 

K{r) oc e-^(2rM)-l^ " 2^ + 0{M-')^ , (A2) 

with Suj'^ = up' — rrijj. Imposing /lo(-R) = gives Slu"^ = 

j. Gauge field dispersion relation The computation is the same as for the Higgs, the 
function ag(r) being the same as /io(r) replacing mn by mw- Imposing a'g{R) = implies 
0^ = + Note that, form equations ([501 HH]) this implies a"(i?) = and ao(i?) = 0, 
which are the correct boundary conditions for these fields. 

k. High energy modes for the Higgs Using the expansion ( lAll) . we see that for u ^ M, 

h^{r) ~ e^™. (A3) 
One of the requirement h{R) = or h\R) = are satisfied if Ru = kn/2, which gives 

APPENDIX B: COUPLINGS IN THE EFFECTIVE ACTION 

In this section, we compute the parameters /. We compute exactly the parameters related 
to the quadratic part of the action, first /„. We start from the normalisation condition for 
ttQ ( 153|) . integrate by part, use the equation of motion (HOi) for ag and the relation 



I dr^ (a>,-" - J(a^)'(r)(a,-")'(r)) 



""'^ I drrAla'^irV^ 
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From which we get 



rfrrA|a^(r)|2 = - -. (Bl) 



This, together with the normahsation condition ( l53l) gives us 



To calculate //i, we integrate by part and use the equation of motion fl33l) 

4 = ^ y rfrrA(^(/i'^)'(r)(/i-7(r) + ;^|/i"(r) 
rfrrA(02 -m^)|/i'^(r) 



^2 



2r2 

^2 _2 \ I i,nl„\\2 



which, using the normalisation condition fl59l) and the dispersion relation, gives J/i = 1. 

The other parameters Ib,Ic,Id,Ie are computed to (9((M_R)~^, (m/^i?)"^) using the 
asymptotic expansion (lAll) for the mode functions. Indeed because of the presence of the 
warp-factor and the fact that the modes , /i" are roughly constant, the integrals Ib,c,d,e 
are dominated by values close to the brane (r = R). At these points, the argument of 
the Laguerre functions is large and the asymptotic expansion can be used. We first com- 
pute Id,Ie dHS]), to this aim, we first have to normalise h^{r). From flA2p . and imposing 
J (irrA|/i"(r)p = 1, we get 



1 M 



^MiR-r)-rn„ ^"^^ ' x i^!" „ 1 „^ (2rl^H) (B3) 




e-^'-"->^ 1 + O 



R \r J V V MR 

Some straightforward algebra gives 

Id = [ dr Ar K^K^h"'-"' = W ^ ( 1 + C ( -— ) ) , (B4) 



R V \MR 



Note that these integrals, at least at leading order, do not depend on the values of n, m, k. 

To compute the remaining /;,, /c, we need the asymptotic expansion for a^{r). It only 
differs from /i"(r) by the normalisation condition. As a^{r) is constant near r = R, we can 
neglect the term Kctg P in the normalisation condition flS51) . which becomes J (ir^|agp = 
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1. This implies that 



„2 



^(f)'"^e-<->^(l + o(-i^)). (B6) 



Using (1B3IIB6P and performing some straightforward algebra gives 

V l\. ( 71 TT) 1 1 



and 



R V XMR'mlrR^ 





rA 




2 




rA 


dr 


2 


H 


+ 




'r 



r2 ^ V Vm2^i?2 



APPENDIX C: COUPLINGS AT VERY HIGH ENERGY 



We already computed the effective coupling A = Ao,o,o,o = MX in the low energy effective 
action fITOl) but we are now interested in couplings XkiMMM^ with ki large. 

At very high energy (0^ ^ M), the wave functions h'^{r) are large near the centre of the 
disk (r = 0). The integrals for the effective couplings will be saturated by values close to 
r = 0, therefore we make use of the small argument expansion of the Laguerre functions 

^^^'^ r(z. + i) VrXATI) f(AT^ + --J- ^^^^ 

This expansion is introduced in the wave function h^: 



1 M 
'2 2(lu 



v^r(2n + l)r (-n - ^ + i 
While computing the effective action, we will identify 



SttR 



A [ vAdrflhli 

i=i 



with 7(0, 0, 0, 0) = 1 and 7(^1, n2, n^, 77.4) of order unity. 
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7(^1, ^2,^3, 714), (C3) 



APPENDIX D: VORTEX SITTING ON THE BRANE, NUMERICAL ANALYSIS 



To simplify slightly the analysis, we notice that the vortex solution has size of a few 
l/rriH which is in our assumptions M ^ itlh ^ 1/R much smaller that the radius of the 
disk. Furthermore what happens away from the brane is damped by the warp factor. We 
can therefore consider that the disk looks like a half plane in the region where the vortex 
lies. The vortex on the brane still possesses one discrete symmetry: the reflection along the 
perpendicular to the brane passing through its centre. We will therefore solve the equations 
on a quarter of plane. To get a better resolution of the instanton centre, we parametrise the 
quarter of plane in polar coordinates r g]0, oo[, 9 G [0, 7r/2]. 

We choose the Aq = gauge and use the remaining time independent gauge freedom to 
cancel the radial component of the vector field'': Ar = 0. We also use dimensionless variables 
as in (fT9|) . The energy functional reads: 



P 
+— 
2 



{drcyf + ( -dea - A 




(Dl) 



Numerically, we minimise the energy functional using as starting configuration, the usual 
vortex configuration A(r), /(r) from fll6lll7p . 

A{r,e)=A{r), f{r,e) = f{r), a{r,9) = 9. (D2) 

The convergence to the exact solution is achieved by introducing an artificial time depen- 
dency and solving Hamiltonian equations: 

d 5E d 6E d _ 6E 

dt'-'JIe' Jt~'6H' Jf-~J^- ^^^^ 

Note that an upper bound on the vortex energy can be found integrating the energy func- 
tional (1D1I) with the usual vortex configuration in a flat 2+1 dimensional space (1161117^ . 

To avoid trivial divergences at r = 0, we start the integration at r = e ~ 10~^ and 
consider a flnite radius R to end the integration.* As we use polar coordinates to describe 
the Higgs fleld, we should also arrange the boundary conditions to avoid that H vanishes 



''' We are dealing with static (time independent) configuration 

^ In dimensionless units we take R — 5. We checked that this do not introduce sizable errors. 
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near r — e and rescale the field a such that it is single valued at r = 0, a{r, 9) — > ^/ra{r^ 9), 
this allows to fix correctly the value of a so that the action remain finite^ at r — > 0. Boundary 
conditions are set on the four edges; 

A0(s,9)^a(s)r^O, H(s,9) ^ f(s) r^O, a(s,9) ^ 9V£ ^ 0, (D4) 
deAe{r,0) = 0, 9ei/(r, 0) = 0, a(r, 0) = 0, (D5) 
A0{R,9) + RdrA0{R,9) ^0, H{R,9)^-, a{r,9) - 2Rdra{R,9) ^ 0, (D6) 

deAe{r,^) = Q, deH{r,^) = 0, a(r,|) = |V^. (D7) 

The conditions at r = and r — R are chosen to get a finite action and the conditions at 
^ = and 9 — 7r/2 are dictated by the discrete symmetries and the topological charge. 
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